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I. INTRODUCTION 



Bottomonium production at high energies has provided a useful laboratory for testing 
the high-energy limit of quantum chromodynamics (QCD) as well as the interplay of per- 
turbative and non-perturbative phenomena in QCD. The factorization formalism of non- 
relativistic QCD (NRQCD) [l| is a rigorous theoretical framework for the description of 
heavy- quarkonium production and decay. The factorization hypothesis of NRQCD assumes 
the separation of the effects of long and short distances in heavy-quarkonium production. 
NRQCD is organized as a perturbative expansion in two small parameters, the strong- 
coupling constant and the relative velocity v of the heavy quarks. 

The phenomenology of strong interactions at high energies exhibits a dominant role of 
gluon interactions in quarkonium production. In the conventional parton model 11^, the 
initial-state gluon dynamics is controlled by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equation P|. In this approach, it is assumed that S > fi^ ^ ^qcd' where 
is the invariant collision energy, fx is the typical energy scale of the hard interaction, and 
Aqcd is the asymptotic scale parameter of QCD. In this way, the DGLAP evolution equation 
takes into account only one big logarithm, namely ln(/i/AQCD)- In fact, the collinear-partron 
approximation is used, and the transverse momenta of the incoming gluons are neglected. 

In the high-energy limit, the contribution from the partonic subprocesses involving t- 
channel gluon exchanges to the total cross section can become dominant. The summation 
of the large logarithms ln(-\/S'/ j^) in the evolution equation can then be more important 
than the one of the ln(yu/AQCD) terms. In this case, the non-collinear gluon dynamics 
is described by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation Q]. In the 
region under consideration, the transverse momenta (fcr) of the incoming gluons and their off- 
shell properties can no longer be neglected, and we deal with reggeized t-channel gluons. The 
theoretical frameworks for this kind of high-energy phenomenology are the /cr-factorization 
approach 0, 0] and the quasi- multi-Regge kinematics (QMRK) approach p, 0], which is 
based on effective quantum field theory implemented with the non-abelian gauge-invariant 
action, as suggested a few years ago [9|. Our previous analysis of charmonium production 



at high-energy colliders using the high-energy factorization scheme jlOl llll | has shown the 
equivalence of the /cr-factorization and the QMRK approaches at leading order (LO) in a^. 
However, the /c5"-factorization approach has well-known principal difficulties at next-to- 
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leading order (NLO). By contrast, the QMRK approach offers a conceptual solution of the 
NLO problems In our LO applications, the QMRK approach yields similar formulas as 
;he fcr-fectoriztion approach, so that we can essentially continue using our previous results 
11 1 obtained in the /cT-factorization formalism using the Collins-Ellis prescription j^. 



IG 



This paper is organized as follows. In Sec. IHl the QMRK approach is briefly reviewed. In 



Sec. IIIH we explain how the analytic results of Refs. [10|, lll[ relevant for our analysis may be 
converted to the QMRK framework. In Sec. II VI we perform fits to the transverse-momentum 
{Pt = IPtI) distributions of inclusive bottomonium production measured at the Fermilab 
Tevatron to obtain numerical values for the non-perturbative matrix elements (NMEs) of 
the NRQCD factorization formalism. In Sec. we summarize our results. 



II. QMRK APPROACH 

In the phenomenology of strong interactions at high energies, it is necessary to de- 
scribe the QCD evolution of the gluon distribution functions of the colliding particles start- 
ing from some scale /io, which controls a non-perturbative regime, to the typical scale n 
of the hard-scattering processes, which is typically of the order of the transverse mass 
Mt = \J M'^ + |prP of the produced particle (or hadron jet) with (invariant) mass M 
and transverse two-momentum py. In the region of very high energies, in the so-called 
Regge limit, the typical ratio x = fi/y/S becomes very small, a; ^ 1. This leads to large 
logarithmic contributions of the type [a^ ln(l/x)]'^ in the resummation procedure, which is 
described by the BFKL evolution equation |^ for an un-integrated gluon distribution func- 
tion $(x, |qTp,/i^), where Qt is the transverse two-momentum of the gluon with respect to 
the flight direction of the incoming hadron from which it stems. Accordingly, in the QMRK 
approach [7] , the initial-state t-channel gluons are considered as reggeons (or reggeized glu- 
ons). They carry finite transverse two-momenta Qt with respect to the hadron beam from 
which they stem and are off mass shell. 

Reggeized gluons interact with quarks and Yang-Mills gluons in a specific way. Recently, 
in Ref. jsl, the Feynman rules for the effective theory based on the non-abelian gauge- 
invariant action [9] were derived for the induced and some important effective vertices. The 
induced vertex for the transition from a reggeized gluon to a Yang-Mills gluon — > g (PR 
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vertex) shown in Fig. CJa) has the form: 



(1) 



where (n~^ 



P2 /E2, Pi 2 are the four-momenta of the colhding protons, 



and £'1,2 are their energies in the center-of-mass frame. We have (n 



±\2 



0, 



n 



2, 



and 5* = (Pi + P2) = AEiE2. For any four-momentum k'^, we define k = k ■ n . It is easy 
to see that the four-momenta of the reggeized gluons can be represented as 



12 



0. 



(2) 



The induced interaction vertex of one reggeized gluon with two Yang-Mills gluons (PPR 
vertex) depicted in Fig. [Tfb) reads 

t 
kt 



(3) 



where Qs = \/47m^ is the gauge coupling of QCD. The reggeized- gluon propagator displayed 
in Fig. m^c) is given by 



(4) 



The Lagrangian of the effective theory [9j also contains the standard gluon-gluon and quark- 
gluon interactions for Yang-Mills gluons. 

Using the Feynman rules for the induced vertices (0) and (jH)) and the ordinary vertices, 
one can construct effective vertices, which obey Bose and gauge symmetries. For example, 
the effective three-vertex that describes the production of a single Yang-Mills gluon with 
four-momentum k^ = qi + and color index h by two-reggeon annihilation + — > g 
(PRR vertex) shown in Fig. |21 reads 

^tbcYili^ k, q2) = Vf^Yi-qi, -q2, k){n+)p{n-)^ + T''~^^{qi, q2, k){n+)p + T''jY{q2, qi, k){n')„ 



^Qsf 



cba 



[n 



V i„^^sr\- i„+v (g- + 41 + (91 - 92)" 



92 



(5) 



where 



V^YYik,, k2, k,) = -g,r'' Uki - k2rg^^ + {k2 - k,fg^^ + {k, - k^Yg 



(6) 



is the Yang-Mills three-gluon vertex, with all four-momenta taken to be outgoing, and we 
exploited the following relation 

The gauge invariance of the effective theory P| leads to the following condition for am- 
phtudes in the QMRK: 



lim \A{R + R ^ n + X)]"^ = 0. (8) 

In the QMRK approach, the hadronic cross section of quarkonium (Ti) production 
through the process 

p + p^n + X (9) 
and the partonic cross section of the two-reggeon fusion subprocess 

R + R-^n + X (10) 

are related as 

J Xi J IT ^ ^ J X2 J 7( 

^(x2^ci2T\\^^'')da{R + R^n + X). (11) 



X 



where $ (x, |qTp,/i^) is the un-integrated gluon distribution function in the proton, Xi = 
qi/{2Ei) and X2 = q2/{2E2) are the fractions of the proton momenta passed on to the 
reggeized gluons, and the factorization scale /i is chosen to be of order Mt- The collinear 
and un-integrated gluon distribution functions are formally related as 

xG{x, /i^) = J^^ c/|qT|'$ (x, M\ /i') , (12) 

so that, for qi^ = q^2T = 0, we recover the conventional factorization formula of the collinear 
parton model, 

da{p + p H +X) = J dxiG{xi, jj^) J dx2G{x2, n'^)da{g + g ^ H + X). (13) 

The partonic cross section or process (fTUI) may be evaluated as 

da{R + R^n + X)= ^ \A{R + R^l-L + X)\'^d^, (14) 

ZX1X20 



where 2x1X28 is the flux factor of the incoming particles, A[R+R — > Ti+X) is the production 
amplitude, the overbar indicates average (summation) over initial-state (final-state) spins 
and colors, is the phase space volume of the outgoing particles, and 



Af^ , '"'"f , (15) 
i6|qiTr|q2T| 



is a normalization factor that ensures the correct transition to the collinear-parton limit. 
This convention implies that the partonic cross section in the QMRK approach is normalized 
approximately to the cross section for on-shell gluons when qiT = q2T = 0. 

In our numerical calculations, we use the un-integrated gluon distribution functions by 
31iimlein (JB) Q|, by Jung and Salam (JS) and by Kimber, Martin, and Ryskin (KMR) 
Ifil ]. A direct comparison between different un-integrated gluon distributions as functions 
of x, Ikj^P, and /i^ may be found in Ref. Note that the JB version is based on the 

BFKL evolution equation 0| . On the contrary, the JS and KMR versions were obtained using 
the more complicated Catani-Ciafaloni-Fiorani-Marchesini (CCFM) evolution equation jisl ]. 
which takes into account both large logarithms of the types ln(l/x) and ln(yu/AQCD)- 



III. RELATION BETWEEN QMRK AND fcT-FACTORIZATION APPROACHES 



In this section, we obtain the squared amplitudes for inclusive quarkonium production 
via the fusion of two reggeized gluons in the framework of QMRK jsl and NRQCD Q]. We 



work at LO in and v and consider the following partonic subprocesses 



R + R 
R + R 



H 

n 



3p(l) 3c(8) lc(8) 3d(8) 



J 5 '-'1 

+ 9- 



^0 ' -'J 



(16) 
(17) 



This formalism also allows for a consistent treatment at NLO, which is, however, beyond 
the scope of this paper and needs a separate investigation. 

According to the prescription of Ref. 1^, the amplitudes of processes fITBj) may be obtained 
from the five Feynman diagrams depicted in Fig. 01 Of course, the last three Feynman dia- 
grams in Fig. 01 can be combined through the effective PRR vertex. The Feynman diagrams 
pertinent to process (fT7|l are shown in Fig. |31 

The LO results for the squared amplitudes of subprocesses (fT^ and (fT7|) that we find 



by using the Feynman rules of Ref. j8[ coincide with those we obtained in Ref. 



Ill ] in the 



6 



/cT-factorization approach. The general relation between the squared amplitudes in both 
approaches is 



N\A{R + R^H + X)\' = \A^^{R + R^n + X)\''. (18) 

The formulas for the 2 — > 1 subprocesses are listed in Eq. (27) of Ref. In the case of 
the 2 — > 2 subprocess (fTTjl . our analytic results were not included in the journal publication 
of Ref. Uj for lack of space. However, they are given in Eq. (38) of the preprint version of 
Ref. and may be obtained in FORTRAN or Mathematica format by electronic mail upon 
request from the authors. 

The differential hadronic cross section of process (fTTjl may then be evaluated from the 
squared matrix elements of processes (|TT)|) and (fTTj) as indicated in Eqs. (46) and (48) of 
Ref. Q. 

IV. BOTTOMONIUM PRODUCTION AT THE TEVATRON 

The CDF Collaboration at the Tevatron measured the pt distributions of T(IS'), T(2S'), 
and T(35') mesons in the central region of rapidity (y), \y\ < 0.4, at ^/S = 1.8 TeV (run I) 
jigj l and that of the T(IS') meson in the rapidity regions \y\ < 0.6, 0.6 < \y\ < 1.2, and 
1.2 < \y\ < 1.8 at = 1.96 TeV (run II) [20]. In both cases, the S-wave bottomonia were 
produced promptly, i.e., directly or through non-forbidden decays of higher- lying S- and P- 
wave bottomonium states, including cascade transitions such as T(3S') X6i(2P) T(IS'). 

As is well known, the cross section of bottomonium production measured at the Tevatron 
at large values of pr is more than one order of magnitude larger than the prediction of the 

linear parton model [22i]. Switching 



color-singlet model (CSM) j21| implemented in the co^ 

rom the CSM to the NRQCD factorization formalism 

23l | somewhat ameliorates the situation in the large-pT region, at pt ^ 10 GeV, but still does 
not lead to agreement at all values of px- On the other hand, the shape of the px distribution 
can be described in the color evaporation model b4^ improved by the resummation of the 
large logarithmic contributions from soft-gluon radiation at all orders in in the region of 
Pt < M ^'2^ . However, the overall normalization of the cross section can not be predicted 



l| within the collinear parton model 



in this approach 



25|. 



In contrast to previous analyses in the collinear parton model, we perform a joint fit 
to the CDF data from run I J^] and run II [2^ for all pr values, including the small- 
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Pt region, to extract the color-octet NMEs of the T{nS) and Xbji'^P) mesons using three 
different un-integrated gluon distribution functions. Our calculations are based on exact 
analytical expressions for the relevant squared amplitudes, obtained in the QMRK approach 
as explained in Sec. Illll 

For the reader's convenience, we list in Table |1] the inclusive branching fractions of the 
feed-down decays of the various bottomonium states, which can be gleaned from Ref. 
Theses values supersede those presented in Ref. 



in Refs. 



Since the T{nS) mesons are identified 



19|,|20| heir decays to pairs, we have to include the corresponding 

branching fractions, which we also adopt from Ref. j^], B(T{1S) fi^ + ji^) = 0.0248, 
B{T{2S) + = 0.0131, and B{T{3S) + = 0.0181. We take the pole mass 

of the bottom quark to be = 4.77 GeV. 

We now present and discuss our numerical results. In Table |nj we list our fit results for 
the relevant color-octet NMEs for three different choices of un-integrated gluon distribution 
function, namely JB Q, JS Q, and KMR Q- 

The relevant color-singlet NMEs are not 
fitted. The color-singlet NMEs of the T{nS) mesons are determined from the measured 
partial decay widths of T{nS) + using the vacuum saturation approximation and 

heavy-quark spin symmetry in the NRQCD factorization formulas and including NLO QCD 
radiative corrections 



28[. The partial decay widths of xm(^-P) 27, from which the color- 
singlet NMEs of the Xbj{nP) mesons could be extracted, are yet unknown. However, these 
NMEs can be estimated using the wave functions evaluated at the origin from potential 
models as was done in Ref. [2^1. We adopt the color-singlet NMEs of the Xbo{nP) 
mesons from Ref. [3]. 

We first study the relative importance of the various color-octet bb Fock states in direct 
T{nS) hadropro duct ion. Previous fits to CDF data were constrained to the large-pT 
region, Pt ~ 8 GeV, and could not separate the contributions proportional to {0'^^^^'>[^Sq^^]) 
and ((9^'-""^^[^Po^^])- Instead, they determined the linear combination 



r 



2 

m 



3 p{8) 



for the value of r that minimized the error on Mj'^"'^\ By contrast, the QMRK ap- 
proach allows us to cover also the small-pT region and thus to fit and 
i/(>^T(n5)|-3p^{8)j^ ggpara^gly^ thanks to the different pt dependences of the respective contribu- 
tions for Pt ~ 8 GeV. This feature is nicely illustrated for direct T(IS') hadroproduction in 



Fig. El where the shapes of the contributions proportional to (O^^^'^^pS'i^'']), {0'^^^^''[^Sq^^]), 
and (O^'-^'^^PPq^'']) are compared. Notice that the peak positions significantly differ, by up 
to 2 GeV. Apparently, this suffices to disentangle the contributions previously combined by 

Eq. (uni). 

In Figs. El 13 and IHl we compare the CDF data on prompt T{nS) hadroproduction in 
run I Il9f with the theoretical results evaluated with the JB 

Q, JS Q, and KMR Q 

un-integrated gluon distribution functions, respectively, and the NMEs listed in Table ITTl In 
each case, the color-singlet and color-octet contributions are also shown separately. Except 
for the JB and KMR analyses of T(3S') production, the color-octet contributions are always 
suppressed, especially at low values of pr- In the JS analysis, the T(IS') and T{2S) data are 
significantly exceeded by the color-singlet contributions for pt ~ 10 GeV, which explains 
the poor quality of the fit, with x^/d.o.f. = 27. In the JB analysis, this only happens for 

^ 2 GeV, so that the value of x^/d-o.f. is lowered by one order of magnitude, being 
X^/d.o.f. = 2.9. By contrast, the KMR gluon yields an excellent fit, with x^/d.o.f. = 0.5, 
and will be the only one considered in the following discussion. Comparing the color-singlet 
and color-octet contributions in Fig. |H1 we observe that the latter is dominant in the T(35') 
case and in the T{2S) case for px ^ 13 GeV, while it is of minor importance in the T(IS') 
case in the whole px range considered. The latter feature is substantiated by the run-II data 
and is refiected in all their y subintervals, as may be see from Fig. ??. 

Notice that the contributions to prompt T{nS) hadroproduction due to the feed-down 
from the Xbj{^P) mesons have been neglected above, simply because the latter have not 
yet been observed and their partial decay widths are unknown. In the remainder of this 
section, we assess the impact of these contributions. For the color-singlet NME, we use the 
potential model result (C»^«(3P)ppa)j^ ^ 2.7 GeV^ H- By analogy to the KMR fit results 
for (CX''o{iP)[35(8)]) and (C»^«(2P)p^(8)]^ TableHH we expect the value of (C^^^^^p)^^^]^ 
to be negligibly small, compatible with zero. Looking at Table HI a naive extrapolation 
from the Xbj{^P) and Xbj{'2P) states suggests that the inclusive branching fractions for the 
Xfej(3P) decays into the T(3S'), T(2S'), and T(IS') states could be about 12%, 9%, and 7%, 
respectively. These decays generate further cascade transitions, whose inclusive branching 
fractions follow from these estimates in combination with the entries of Table H] Including 
all these ingredients, we repeat our KMR fit to the CDF data. As illustrated in Fig. [TUl 
for prompt T{nS) hadroproduction in run I, the CDF data can be fairly well described 



9 



in the QMRK approach to the CSM, while the color-octet contributions turn out to be 
negligibly small. We note in passing that a similar observation, although with lower degree 
of agreement between data and theory, can be made for the JB gluon, while the JS gluon 
badly fails for px ~ 10 GeV. 



V. CONCLUSION 



Working at LO in the QMRK approach to NRQCD, we analytically evaluated the squared 
amplitudes of prompt bottomonium production in two-reggeon collisions. We extracted 
the relevant color-octet NMEs, (C^psf^]), {0'^[^Sf^]), and {0'^[^P^^^]) for H = T{1S), 
T{2S), T{3S), Xboi^P), and xw(2-P), through fits to pt distributions of prompt T{nS) 
hadroproduction measured by the CDF Collaboration at the Tevatron in pp collisions with 
75 = 1.8 TeV Q and 1.96 TeV Q using three different un-integrated gluon distribution 
functions of the proton, namely JB Q, JS H, and KMR Q. The fits based on the 
KMR, JB, and JS gluons turned out to be excellent, fair, and poor, respectively. They 
yielded small to vanishing values for the color-octet NMEs, especially when the estimated 
feed-down contributions from the as-yet unobserved X6j(3-P) states were included. 

The present analysis, together with a recent investigation of charmonium production 
at high energies |llj| , suggest that the color-octet NMEs of bottomonium are more strongly 
suppressed than those of charmonium as expected from the velocity scaling rules of NRQCD. 
We illustrated that the QMRK approach 0, 0] provides a useful laboratory to describe the 
phenomenology of high-energy processes in the Regge limit of QCD. 

LO predictions in both the collinear parton model and the QMRK framework suffer from 
sizeable theoretical uncertainties, which are largely due to unphysical-scale dependences. 
Substantial improvement can only be achieved by performing full NLO analyses. While the 
stage for the NLO NRQCD treatment of 2 — 2 processes has been set in the collinear parton 
model js^, conceptual issues still remain to be elaborated in the QMRK approach. Since, 
at NLO, incoming partons can gain a finite krp kick through the perturbative emission of 
partons, one expects that essential features produced by the QMRK approach at LO will 
thus automatically show up at NLO in the collinear parton model. 
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TABLE I: Inclusive branching fractions of the feed-down decays of the various bottomonium states. 



In\Out 


T{3S) 


Xb2{2P) 


Xbi{2P) 


Xbo{2P) 


T(25) 


Xb2{lP) 


XbiilP) 


XboilP) 


r{is) 


T{3S) 


1 


0.114 


0.113 


0.054 


0.106 


0.00721 


0.00742 


0.00403 


0.102 


Xb2{2P) 




1 






0.162 


0.0110 


0.0113 


0.00616 


0.130 


Xbii^P) 






1 




0.21 


0.0143 


0.0147 


0.00798 


0.161 


X6o(2P) 








1 


0.046 


0.00313 


0.00322 


0.00175 


0.0167 


T{2S) 










1 


0.068 


0.07 


0.038 


0.320 


Xb2{lP) 












1 






0.22 
















1 




0.35 


XboilP) 
















1 


0.06 



T(15) ••• ••• ••• 1 



13 



c, ± 



a, ± 6,1/ 



a, /i 




q a, fi 



b, V 



(b) 



6, V 



FIG. 1: Feynman diagrams pertinent to the (a) PR vertex, (b) PPR vertex, and (c) reggeized-gluon 
propagator given in Eqs. ©, and @, respectively. 




FIG. 2: Feynman diagrams pertinent to the effective PRR vertex given in Eq. 
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TABLE II: NMEs of the T(15), T(25), T{3S), Xboi^P), and Xbo{2P) mesons from fits to CDF 
data from run I 19] and run II in the collinear parton model (PM) 2^ using the CTEQ5L 



27 1 parton distribution functions of the proton and in the QMRK approach using the JB [145, 
15 1, and KMR [l^ un-integrated gluon distribution functions of the proton. The errors on the fit 
results are determined by varying in turn each NME up and down about its central value until the 
value of is increased by unity keeping all other NMEs fixed at their central values. 



NME 


PM [23] 




Fit JB 




Fit JS 




Fit KMR 


(0^(i^)[35f)])/GeV3 


10.9 ± 1.6 




10.9 ±1.6 




10.9 ± 1.6 




10.9 ± 1.6 


^0T{15)[3^(8)]^/QgV3 


(2.0 ±4.1;°:^) X 


10-2 


(5.3 ±0.5) X 10- 


-3 


(0.0 ± 1.8) X 10- 


-4 


(0.0 ±3.1) X 10 


])/GeV3 






(0.0 ±4.7) X 10" 


-4 


(0.0 ±5.2) X 10" 


-5 


(0.0 ±4.3) X 10 


^0T(lS)[3p{8)]^/Qgy5 






(0.0 ± 1.3) X 10" 


-3 


(0.0 ± 1.6) X 10" 


-4 


(9.5 ±2.0) X 10 




(l.4±0.7+o;°) X 


10-1 


(0.0 ±7.6) X 10- 


-4 


(0.0 ±8.7) X 10- 


-5 


(2.1 ±0.9) X 10 




2.4 ±0.4 




2.4 ±0.4 




2.4 ±0.4 




2.4 ±0.4 


(OXi,o(iP)[3^(8)])/GeV3 


(l-Sil-l+Lo) X 


10-2 


(0.0 ±2.1) X 10" 


-3 


(0.0 ±8.4) X 10- 


-5 


(0.0 ± 1.4) X 10 


^0T{25)[35(l)])/QgV3 


4.5 ±0.7 




4.5 ±0.7 




4.5 ±0.7 




4.5 ±0.7 




(l.6±0.6l°:^) X 


10-1 


(0.0 ±5.9) X 10" 


-3 


(0.0 ±4.1) X 10- 


^4 


(3.3 ±0.8) X 10 


^0T{25)[l5(8)]^/Qgy3 






(0.0 ±9.2) X 10~ 


-4 


(0.0 ±8.3) X 10- 


-5 


(0.0 ±3.7) X 10 


^0T(2S)[3pW]^/GgV5 






(0.0 ±2.6) X 10- 


-3 


(0.0 ±2.8) X 10- 


-4 


(0.0 ± 1.6) X 10 




;-i.i± 1.0+°:^) X 10-1 


(0.0 ± 1.5) X 10" 


-3 


(0.0 ± 1.4) X 10- 


-4 


(0.0 ± 7.2) X 10 


(0Xi,o(2P)[3pa)]^/QgY5 


2.6 ±0.5 




2.6 ±0.5 




2.6 ±0.5 




2.6 ±0.5 


(OX.o{2P)[3^(8)])/GeV3 


{0.8±l.ltl-l) X 


10-2 


(1.1 ±0.4) X 10" 


'2 


(0.0 ±2.8) X 10- 


-4 


(0.0 ± 5.7) X 10 


^0T(35)[35(l)])/QgV3 


4.3 ±0.9 




4.3 ±0.9 




4.3 ±0.9 




4.3 ±0.9 


^0T(35)p5(8)]^/Qgy3 


(3.6 ± 1.9+1;^) X 


10-2 


(1.4 ±0.3) X 10" 


-2 


(5.9 ±4.2) X 10- 


-3 


(1.1 ±0.4) X 10 


])/GeV3 






(0.0 ±2.6) X 10- 


'3 


(0.0 ±8.1) X 10- 


-4 


(0.0 ± 2.7) X 10 


^0T(3S)[3pW]^/GgV5 






(2.4 ±0.8) X 10" 


-2 


(3.4 ±4.2) X 10- 


-3 


(5.2 ± 1.1) X 10 




(5.4±4.3+|i) X 


10-2 


(5.2 ±4.4) X 10" 


-3 


(7.4 ± 10.2) X 10 


-4 


(1.1 ±0.5) X 10 


(C)X50(3P)[3pa)]^/Qgy5 


2.7 ±0.7 




2.7 ±0.7 




2.7 ±0.7 




2.7 ±0.7 


xVd.o.f. 






2.9 




27 




0.5 



15 




R, q2, — R, 92, — R, 92, — 



FIG. 3: Feynman diagrams pertinent to processes (|T6|) . 




R, 92, - g, h R, 92, - g, h 



FIG. 4: Feynman diagrams pertinent to process (|17j) . 
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2 4 6 8 10 12 14 16 18 20 

FIG. 5: Contributions to the px distribution of direct T(IS') hadroproduction in pp scattering with 
\/5 = 1.8 TeV and \y\ < 0.4 from the relevant color-octet states. All distributions are normalized 
to unity at their peaks. 
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(a) 




FIG. 6: pt distributions of prompt (a) T{1S), (b) T(25), and (c) T{3S) hadroproduction in pp 
scattering with -v/S" = 1.8 TeV and \y\ < 0.4 including the respective decay branching fractions 
B{T(nS) — > The color-octet (curve 1) and color-singlet (curve 2) contributions, evaluated 

with the JB j^l un-integrated gluon distribution function, and their sum (curve 3) are compared 
with the CDF data from run I |19l |. 
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FIG. 9: pt distributions of prompt T(15) hadroproduction in pp scattering with y S = 1.96 TeV 
and (a) \y\ < 0.6, (b) 0.6 < \y\ < 1.2, (c) 1.2 < \y\ < 1.8, and (d) \y\ < 1.8 including the decay 
branching fractions B{T{1S) fi^ + fi^). The color-octet (curve 1) and color-singlet (curve 2) 



contributions, evaluated with the KMR |.16| un-integrated g. 
sum (curve 3) are compared with the CDF data from run II 



uon distribution function, and their 
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V 
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(b) 




(c) 



FIG. 10: pt distributions of prompt (a) T(15), (b) T{2S), and (c) T{3S) hadroproduction in pp 
scattering with -v/S" = 1.8 TeV and \y\ < 0.4 including the respective decay branching fractions 
B{T(nS) — > /U^ + n^)- The color-singlet contribution including the estimated feed-down contri- 
butions due to the meson, evaluated with the KMR 16] un-integrated gluon distribution 
function, is compared with the CDF data from run I |1^]. 
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